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^^ ' The possibility of treating colour in one-loop amplitude calculations alike the other quantum 

C^^ \ numbers is briefly discussed for semi- numerical algorithms based on generalized unitarity 

^~^ . and parametric integration techniques. Numerical results are presented for the calculation of 

^_^ ' virtual corrections in multi-gluon scattering. 

' \ 1 An algorithmic solution based on generalized unitarity for the evaluation of 

^' colour-dressed one-loop amplitudes 

^: 

'""'■ The evaluation of QCD higher-order corrections in particular for multi-particle processes is in- 

,—1 . dispensable for a good understanding of total and differential cross sections at hadron colliders 

^ \ such as the Tevatron or the LHC. The automation of these calculations at the next-to-leading 

order (NLO) in the strong coupling constant has been a very appealing goal ever since the 
Q^ . Monte Carlo tools dealing with the computation of tree-level cross sections have matured and 

^N \ been greatly optimized. Over the past few years a number of subtraction codes has been de- 

^O ' veloped handling the cancellations of singularities occurring in NLO calculations in a general 

way. Accordingly a standard interf ace has lately been agreed on between Monte Carlo tools and 
one-loop matrix-element programs'^'J^. Recently these programs have also shown a remarkable 
advancement. New techniques based on combining generalized unitarity and parametric integra- 
tion methods have become available and made the computation of multi-leg one-loop corrections 
feasible. Prominent examples are given by the vector boson plus 3 jet NLO predictions provided 
_C^ \ by the BlackHat and Rocket groupg "^ ! 1^1"! and the ti plus 1 jet NLO result^^. These calcula- 

tions separate the treatment of the colour quantum numbers from the other degrees of freedom, 
which makes it harder to fully automatize these approaches. In a recent publicationl^ it has been 
shown that the Ellis-Ciele-Kunszt-Melnikov algorithm'^ can be extended to treat colour along 
the same line with the other quantum numbers. Given the implementation for ordered gluon 
one-loop amplitudes'^, the important changes are: the sums over ordered cuts present in the 
decomposition of the one-loop integrands and amplitudes are replaced by sums over unordered 
partitions including their non-cyclic and non-reflective permutations. Gluon bubble contribu- 
tions now come with a symmetry factor of 1/2!. The integrand's residues are calculated through 
products of colour-dressed tree-level amplitudes obtained from dressed recursion relations where 
internal colour degrees of freedom are summed over. Finally, to extract the lower-point coef- 
ficients, the subtraction terms due to higher-cut contributions are identified by a de-pinching 
procedure that may involve loop-momentum shifts by external momenta. 
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Figure 1: Double-, single-pole and finite-part accuracy distributions obtained from double-precision computations 

of one-loop amplitudes for A'^ = 4, 5, 6, 7 gluons with polarizations A^ — H \ ... and randomly chosen 

non-vanishing colour configurations. Results from the colour-dressed algorithm are compared to those of the 

colour-ordered method (labeled "ordr") indicated by dashed curves. The mean accuracies and the number of 

randomly picked phase-space points (subject to the cuts detailed in the text) are displayed at the top and bottom 

left of the plots, respectively. Phase-space points required to be calculated at higher precision were vetoed. 

2 Numerical calculation of multi-gluon one-loop corrections 



Taking the example of multi-gluon scattering, the performance and major results of the algorithm 
briefly introduced in Section [T] are discussed below. A more comprehensive description of the 
outcomes of these numerical calculations can be found in the original paperl^. 

The convenient way to cross-check the results of the colour-dressed algorithm for the full 
one-loop amplitudes for N gluons is to make use of the colour decomposition approach. Given 
the polarization states A^ G {-|-, — } of the k = 1, . . . ,N gluons and their colours {ij)k using 
the colour-flow notation ik,jk G {1)2,3}, one sums up all relevant ordered amplitudes after 
having multiplied them by their corresponding colour factors. The colour-ordered A^-gluon one- 
loop amplitudes are computed with the algorithm^LiiJ based on the original Ellis-Giele-Kunszt- 
Melnikov method. The quality of the one-loop amplitude determination can be analyzed by 
means of the logarithmic relative deviations of the double (dp) and single poles (sp) and the 
finite part (fp). They are defined as follows: 
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Two independent solutions denoted by [1] and [2] are used to test the accuracies of the single 
poles and finite parts. The analytic structure of the double poles is particularly simple and pro- 
portional to the Born amplitudes, A^jp ^^ = — (cr/e^) N Nq M.^^' where A^'c denotes the number 
of colours. This allows for a direct comparison with the semi-numerical results. Figure [T] shows 
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Figure 2: Times tat to compute two A'^-gluon one-loop amplitudes for a hard, a simple and random colour 
configurations using 3.0 and 2.66 GHz Intel Core2 Duo processors. The solid and dashed curves represent fit 
results for the colour-dressed (drss) and colour-ordered (ordr) approach, respectively. The times to compute the 
cut-constructible part of the amplitudes using the 4-dimensional (4D) versions of the algorithms are also shown. 



the e-distributions in absolute normalization obtained by using double-precision computations 
for various numbers of external gluons with alternating polarization states. The phase-space 
points are accepted if the generated momenta satisfy the cuts: \r\n\ < 2 , p_|_^„ > 0.1 \Ei + £"2! 
and Ai?fc„ > 0.7 where ry„ and p±^n respectively denote the pseudo-rapidity and transverse mo- 
mentum of the n-th outgoing gluon (n = 3, . . . , N); the l^Rkn describe the pairwise geometric 
separations in rj and azimuthal-angle space of gluons k and n. For N > 5 gluons, a small fraction 
of events requires higher than double precision calculations to reliably determine the master- 
integral coefficientsO These phase-space points can be identified by a simple procedure testing 
the stability of the solutions for the bubble coefficients that contribute to the cut-constructible 
part. As seen in the plots of Figure [1] the gluon one-loop amplitudes can be determined quite 
accurately for the bulk of the events. The tails are sufficiently under control and fall off more 
steeply for larger e-values owing to the veto on points that yield unstable solutions in double- 
precision calculations. The distributions and peak positions of the double poles are rather stable 
while those of the single poles and finite parts noticeably shift to larger e with an increasing 
number N of gluons. In all cases the dressed approach is seen to provide more accurate results 
than the method relying on the colour decomposition. 

The complexity of the colour-dressed tree-level recursion relations and the asymptotic be- 
haviour of the Stirling numbers of the second kind that govern the growth of the number of the 
unitarity cuts both follow an exponential scaling law when more legs are added. The colour- 
dressed algorithm is hence expected to scale similarly with A^, i.e. the computing time tn is 
proportional to x where x is an attribute of the implemented algorithm. Relying on the colour 
decomposition, the ordered one-loop amplitude calculations exhibit polynomial complexity while 
the number of orderings that need be evaluated for one colour configuration grows factorially. 
The effective growth however will be determined by the increase in the number of non- vanishing 
orderings. Figure [2] depicts how the A^-gluon one-loop amplitude computation times scale with 
increasing A^ for different examples of assigning the gluons' colour states. Amplitudes with a 
small number of gluons, A^ < 6, can always be calculated faster using the colour decomposition. 
This remains true as long as there is a sufficient number of vanishing orderings, i.e. one deals 
with simple colour configurations only. In all other cases - for random colour assignments most 
importantly - the colour-dressed method takes over for N > 7 owing to its milder exponential 
growth. The curve fitting for the general case indicates base-values of x = 7.3(1) for the dressed 
algorithm versus x = 9.5(1) when using the colour decomposition. One notices that the factorial 
growth caused by the sum over orderings has been tamed to become effectively exponential. 



For A'^ = 6 and 7 gluons, about 3% and 10% of the events need be treated in higher precision, respectively. 
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Figure 3: Consistency check for and relative errors of the phase-space integrations of colour-summed and colour- 
sampled amplitudes at the Born plus virtual-correction level. The distributions are given as functions of the 
number of generated uniform phase-space points; cuts are specified in the text. After 15900 events one obtains 
('S'mc^'>/(S'^oi^^^) = 0.939 ±0.039(MC)±0.028(col) with the left plot showing errors due to colour sampling only. 



Utilizing the ability of calculating the full one- loop amplitudes for multiple gluons, one can 
study the effect of the Monte Carlo sampling over non- vanishing ("Non-Zero") colour configu- 
rations on the performance of phase-space integrations for the Born plus virtual contributions. 
Figure [3] shows the simplest example of 2 — > 2 gluon scattering. The colour-sampled and colour- 
summed integrals are defined as S\j^q = Wcoi x ^ and 5^^^ = X^coi ^) respectively, where the 
colour weight W^oX depends on the actual colour configuration ( "col" ) and the kernel is given by 
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with Qo = 0.12 



(2) 



In the left panel of Figure [3] the sampled and summed integrations are shown to converge after 
about 10^ Monte Carlo steps, NmCi while in the right panel their relative errors are displayed as 
a function of Nmc- One finds that the colour sampling only introduces an additional integration 
uncertainty, which can be reduced more easily than the full colour sum can be carried out. 
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